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Abstract 

We show that many integrals containing products of confluent hypergeometric 
functions follow directly from one single integral that has a very simple formula in 
terms of Appell's double series F 2 . We present some techniques for computing such 
series. Applications requiring the matrix elements of singular potentials and the per- 
turbed Kratzer potential are presented. 
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1. Introduction 

Landau and Lifshitz [1] have discussed an important type of integral containing a 
pair of confluent hypergeometric functions 1F1, namely 

oc 

j; p (a,a')= J 'f- 1+s e- ht 1 F 1 {a;r,kt)iFi{a';'y-p;k't)dt, (1.1) 
o 

where s and p are positive integers, and values of the parameters are supposed such 
that the integral converges absolutely; Re(h) > . These authors mentioned that a 
general formula for such integrals can be derived by means of a method proposed 
earlier by Gordon [2], but this is so complex that it cannot be conveniently used. 
Instead, they investigate the particular case J°°(a, a') and they show that, for h = 

k+k' 
2 > 

Ahk 1 

J 7 00 (a, a') = 2^r( 7 )(fc + k') a+a ~ 7 (fc' - k)- a (k - k'y a 2*1 (a, a'; 7 ; - _ ). (1.2) 

The functions 1.F1 and 2 Fi just mentioned are particular cases of the generalized 
hypergeometric function 

p 

oo n {ai) k xk 

p F g (ai,a 2 , . . .,a p ;(3i,p 2 , ■ ■ ■ ,(3 q ;x) = ^ ^ — , (1.3) 

3=1 

where p and q are non-negative integers, and none of the parameters (j = 
1,2, ...,q) is equal to zero or to a negative integer. It is known that if the series 
does not terminate (that is to say, if none of the a, , i = 1,2, . . . ,p , is a negative 
integer), then, in the case p = q + 1 , the series converges or diverges according as 
\x\ < 1 or \x\ > 1 . For x = 1 on the other hand, the series is convergent, provided 

q p 

Pj - a i > 0- Note (6) n , the shifted factorial (or Pochhammer symbol), is defined 

3=1 i=l 

by 

T(S + m) fi, if m = . . 

(6)m ~ T(S) \5(5 + l)...(5 + m-l), ifm = l,2,... U ' 4j 

Karule [3], in order to overcome the inconvenient approach suggested by Landau 
and Lifshitz, proposed a method of derivation which allowed him to obtain a general 
formula for the integrals (1.1). He showed that 

k, ,„ fe\ ' ^ l(-s - p) m (oi') m h 



J^( a ,a') = (l--)- a {l--)- a 'Y J 



m=0 



1 - — 



(c - p) m ml V fc ' 

(1.5) 

Ef— s — m) r (a) r ( h\ r , kk! 

" (c) r rl l 1 -^ 2 ^ a + r ^ + ro;c + r; ^-^- fc) ) 



r=0 



It should be clear however that the integral (1.1) and consequently (1.5) is convergent 
for |fc| + \k'\ < \h\ . For h = ^f- , both (1.2) and (1.5) should be interpreted in terms of 
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analytic continuations whenever the integral (1.1) is absolutely converge. For example, 
if a and a' are nonpositive integers and since the confluent hypergeometric functions 
i-Fi becomes polynomials, the condition \h\ > \k\ + \k'\ is no longer necessary and the 
condition Re(h) > is sufficient for the convergent of the integral to take a place. 

In the present article, we show that the integral (1.1) is indeed a special case of 
a more general type which arises in the different fields in physics. For example, the 
computation of the radial matrix elements of the electromagnetic interaction between 
the states of a relativistic electron in the Coulomb field of a point nucleus [4], namely 



/■ 



-V^Fi(a;&;fci)iFi(a';&';fc'i) dt. (1.6) 



The matrix element integral [5] for bremsstrahlung or pair production without Born 
approximation is usually expressed in terms of confluent hypergeometric functions 
(1.6). 

For convenience, in section 2, we present an elementary method to find a general 
formula for the integral (1.6) in terms of the Appell's series F 2 (d;a,a';6,6';^,^-) . The 
series that is absolutely convergent if \k\ + \k'\ < \h\ . Furthermore, we show that many 
integrals [6] involving confluent hypergeometric functions, or a product of confluent 
hypergeometric functions with weight measure {i d ~ 1 e~ /lt }, follow immediately from 
the integral (1.6). In section 3, we prove that the result of Karule, (1.5), follows for 
special values of the parameters in the Appell's series F 2 . This allows us to develop 
some new results useful for the computation of Appell's series F 2 . For example, we 
shall give a generalization of the known identity 

F 2 {d;a,a';d,d;x,y) = (l-x)- a (l-y)- a ' 2 F 1 (a,a';d- 



(i- x )(i- y y 



in terms of the Gauss hypergeometric series 2 F X . We also present some techniques 
for computing the Appell's hypergeometric series F 2 , whereby we discuss several 
recurrence relations for F 2 . In section 4, we present some applications involving the 
computation of the matrix elements for singular potentials and for the perturbed 
Kratzer potential. 

2. An integral of general type and some special cases 

Appell's hypergeometric [7]- [10] series F 2 , which is a function of two complex ar- 
guments x, and y , and of five complex parameters a, ai,a 2 ,h, and b 2 (h,b 2 ^ 
0, -1, -2, . . . ), is defined by 

F2(a;a u arM,b 2 ;x,y)=^}^ {bl)Mn m! n! x V " ^ 

where the Pochhammer symbols (S) n are defined by (1.4). The series converges ab- 
solutely for \x\ + \y\ < 1 and, in general, diverges for \x\ + \y\ > 1 . In the present work, 
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we will assume, unless otherwise stated, that 01,02,61, and b 2 are real with 61 and 
b 2 strictly positive. The following lemma can be seen as an integral representation of 
F 2 , as well as a Laplace transform of the product of two confluent hypergeometric 
functions [11]. 

Lemma 1: 

For Re(d) > , and \k\ + \k'\ < \h\ , 



OC 

/ 



t d -V^Fi(o;6;H)iFi(a';6';/c'i) dt = h~ d T(d) F 2 (d;a,a';b,b';^,j) (2.2) 



Proof: From the series representation (1.3) of the confluent hypergeometric series 
i-Fi , we have 



°° °° r OO OO 

/ t d - 1 e- ht 1 F 1 (a;b;kt) 1 F 1 (a';b';k't)dt= / t^e^ ^ ^ 

J J L ™ — n m — n 



(q) m (q / ) n k m k f ^ n+m 



m -^o m!n! 



(J.) (a)m(aOrt fc m fc / /" d+n+m _ 1 ht 

\H (b)m(b% m\n\ J 

Wu-d-rr^ (d)m+n(a)m(a')n fk'\ n fk\ m 

{ ) hh ( b ^% ml n\ \h) [h] (2 - 3) 

where (1) is justify by the condition \k\ + \k'\ < \h\ and the fact that for large x the 
confluent hypergeometric function iFi(o, c, a;) is asymptotic to [12] 

l F 1 (a,c,x) = ^e x x a - c [l + 0(\x\- 1 ] ] 

(2) follows by means of the definition of Gamma function with Re{d) > and Re(h) > 
0, and the use of Pochhammer's identity 

r(d + m + n) = (d + m) n T(d + m) = (d + m) n {d) m T{d). 

The right hand side of (2.3) is Appell's hypergeometric series F 2 which observation 
completes the proof of the lemma. □ 

Unfortunately, the direct use of the series representation (2.1) for the Appell's hy- 
pergeometric series F 2 does not converge for the variables needed in the computation 
of Landau and Lifshitz's integral (1.1), namely h < k + k' . However, there are several 
analytic continuations of Appell's hypergeometric series F 2 available in the literature 
[4], [12] -[14]. In the next section we shall give several new analytic continuation of F 2 
that can be used directly in the computation of the integral (1.1). First, we show that 
lemma 1 can be used directly to prove many standard and non-standard results. The 
Appell's F 2 (d;a,a';c,c';x,y) function reduces to a Gauss hypergeometric function 2 F X 
if a or a' or just one of the variables (x,y) is zero. 
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Lemma 2: 

1. F 2 (d;0, a';c,c';x,y) = 2 F 1 (d, a'; c' , y) . 

2. F 2 (d; a, 0; c, c'; x, y) = 2 Fi (d, a; c, x) . 

3. F 2 (d;a, a';c, c';0, y) = 2 Fi(d,a';c',y) . 

4. F 2 (d; a, a'; c, c'; x, 0) = 2 Fi(d, a; c, x) . 

These allow us to deduce the following known formula, for Re(d) > and \k\ < \h\ , 

oo 

J t^e-^F^a-^kt) dt = h- d T(d) 2 F 1 (d,a;b;^-). (2.4) 
o 

Many identities follow directly from (2.4) by means of the expansion for 2 F\ in terms 
of elementary functions [12], a few of which we collect in Appendix I for use in further 
applications. Many other formulae can be deduced directly from (2.2), for example, 
formulas (19) to (21) in Appendix I. 

Lemma 3: 

1. F 2 (d;a,a';a,a';x,y) = (l — x — y)~ d 

2. F 2 (d; a, a'; d, a'; x, y) = (1 - y) a - d {l -x- y)~ a 

3. F 2 (d; a, a'; a, d; x, y) = (1 - x) a '- d (l -x- y)~ a ' 

Consequently, the following formulae follow: for Re{d) > and |fc| + \k'\ < \h\ , we have 

oo 

J t^e-^F^aiaikthFiia'ia'ik't) dt = T(d)(h - k - k')- 1 (2.5) 
o 

which is expected since iFi(/3;/3;x) = e x . Also, for Re(d) > and \k\ + \k'\ < \h\ , we 
have 

oc 

J f , - 1 e- h *if , i(a;d;fct)iFi(a / ;a / ;fc / t) dt = T{d){h - k') a - d {h - k - k 1 )- 1 . (2.6) 
o 

Recently, Cuyt et al [22] developed a finite sum representation of Appell's hyper- 
geometric series [7]- [10] 

Fl (a; b, 6'; c; *, y) = f) £ Wp^U^ ^,^ (2 . ?) 

m— n— 

under the conditions that the parameters a, 6, 6', and c are positive integers; namely 

Lemma 4: For any non-negative integers a, s, t, d and x ^ y , we have for \x\ < 1 , 
I»I<1, 

r(o + i)r(d + i) 



T(a + d + 2) 



Fi(a + 1, s + l,i + 1; a + d + 2;a;,y) = 



d / , n t 



m=0 v 7 k j=0 v / \U ) 

t-a-m \- f k + t\ {-yf \ 

M < / («-i/)*+ fc + i / 



fc=0 
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where 



a+m 



E 



k=0,k^t-j 
a+m 



T)H)' |i ^ l -(.-: i|( -"""- ti - A 



c- f (°r) '- 1 ' j ll "o 1 + " t -T' - ('-") '-'r^d-x). 

j—O^j^s k 

Here, (j^^j is the binomial coefficient and all sums are understood to be zero where 
they are not defined. 

These results can be used along with the following reduction formulae for F 2 to 
obtain many new identities for the integral (2.1) with positive integer parameters. 

Lemma 5: 

1. F 2 (d;a 1 a';d,c';x,y) = (1 - x)~ a F 1 {a', a, d - a;c'; +^,y) 

2. F 2 (d;a,a';c,d;x,y) = (1 - y)- a ' F 1 (a, d - a', a'; c; x, y+) 

For example; since 



F 2 (4; 2, 1; 4, 2; x, y) = (1 - ^"^(l, 2, 2; 2; 

1 

+ 



2/ 



1-x 
1 



a+l — y) x 2 (l — x)(l — .x — y) x 3 y 



ln(l - y) - ln(l 



1 -x' 



(2.8) 



we have 



oo 

J t z e~ h \F x (2; 4; fci)^ (1; 2;fc'i)di = T(4) | - 



k 3 k' 



) k 2 (h-k)(h-k-k') 



ln(/i - k') + ln(/i - fc) - ln(/i) - ln{h - k - k') 



(2.9) 



3. Calculations of Appell's hypergeometric function F 2 

In this section, we present our main results concerning the computation of F 2 
and its analytic continuation. The integral formulae in this section are mostly new, as 
distinct from the many known results which we can recover, some of which we have 
listed in the appendix. 

Lemma 6: For s and d positive integers, 
1. In terms of Fi , we have 



F 2 {c + s;a,a';c,c-p;^,^-) = ( 1- ^- j 



k'\ a (a') m (-s -p) 7 



m=0 



{c-p) m ml 



k' 



k k 

x F 1 (a;c + s - a' ,m + a';c; -; —). 



(3.1) 
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2. In terms of 2 Fi , we have 

W / u\~ a / u'\ S +P 



„, , , k k' ( k\~ a ( k'Y a ^ \ {a') m {-8-p) m {, h\ 

E fm (a)r(— s — m) r f h\ T ' , . kk' 
, n \c) r rl {'--k) ^(a + r;a> + m;c + r; ( , _ fc , )(fe _ fc) ) 



X 

r=0 



(3.2) 

where \k\ + \k\' < \h\ . 

Proof: From the integral representation [16] of the Appell's hypergeometric function 
F 2 , we have 

k k' T(c)T(c-p) 
F 2 {c+s;a,a';c,c-p;-,-)- 



h' h' T(a)T(a')T(c - a)T(c - p - a') 
l l ^ 

u°- V'-^l - u) c - a ~\l - v) c - p - a '- 1 (l - - 'Lv)- C - S dudv 







under the conditions: c>a>0,c-p>a'>0. We shall look first at the integral 

i i 

J h h h J h — ku 

o o 

Since the integral on the right-hand side is now the integral representation [17] of the 
Gauss hypergeometric function 2 Fi , we have 



k' 

x 2^1 ( c + s, a'; c-p; — )■ 



From Pfaff's transformation [18] of 2-F1 , namely 



2 Fi(a,&;c;;j) = (1 - ^"^(a, c - 6; c; — ^-), (3.3) 

2 — 1 

we have 

s+p 



h-ku' ^ {c-p) m m\ [ -h + ku + k'' 



m=0 



since 2 F 1 on the right-hand side of (1) is a polynomial of order s + p. Therefore, we 
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may now write 



7p / | / k k' T(c) , ^ (a') m (-s -p) m i ,n 

F 2 (c + s;a,a;c,c-p;-,-) = — — — -h + E (r _ n ,^ (-* ) 



o 



'/i' /i r(a)r(c-a) ^-^ (c-p) m m! 
1 

x y w a - 1 (I-u) c - a - 1 (/ l -fcu)- c - ;s+a >-fc U -fc') _m_a 'dw 
o 

r ( C ) / ft W \p («')mH-p)m ( fc / 

r(a)r(c-a)^/i- k'' ^ (c-p) m m\ [ k> - h> 
i 

f u a -\l - uf- 1 - 1 ^ - ^ u )- c - s+a '(l - J^-)- m - a 'du 
J h h — k' 

o 

The last integral is an integral representation [19] of Appell's series Fx , namely 

f , ( « ; ^,^4i)= f ^/,--- (1 -,r-'(4)--'(i-A r 

o 

Therefore 

m— 

This proof (3.1) in the lemma. From the identity [20] 

Fi(a; 6, 6'; c; x, y) = (1 - z)""^ (a; c - 6 - 6', 6'; c; ^-p |-|) 

we have 

l ■ I k k k . _ . k . _ / ^— \ (a)m( s f)m /-, ft \— m 

F 2 (c + W ;c,c- P ;-,-) = (l--) (I--) E (c - p)mm! ^ ~ ^ X 

m— v y 

F l( a; - m, a' + m; c; ^A_; {h _ ™ {h _ k) ) 
As the second parameter of the Fi function is a negative integer, we may now write 



m=0 

s+m 



E(a) r (— s — m) r k r . kk' 

r (c) r {k-hyr\ 2Fl(a + r; ° + m; C + r; (h-kOih-k^ 



because of the identity 



Fx (a; 6, 6'; c; x, y) = ^ ^^f^^ F (a + r,b':c + r; y) 
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which complete the proof of the lemma. □ 

We may also apply Pfaff's transformation (3.3) to the hypergeometric function 
2-Fi in (3.2) which allow us to write (3.2) as follows 

k fc\ ^ fc s _„ /a fc\_„' l(-s -p)m(a') m . h 



F 2 {c + s;a, a'; c, c - p; ^) = (1 - - ^ £ 



(c - p) m ml k 



(I--)-™x 



E Vm (a)J-8-m) r k r , kk' . „/ „, , fcfc' 

o (cMfc-^H ^ - mom)) 



As an immediate result of (3.2) the following identity follows: 

F 2 (c;a,a';c,c;x,y) = (1 - x)~ a (l - y)~ a ' 2 F 1 (a, a'; c; ^ -). (3.4) 

A general formula for Landau and Lifshitz's integral type (1.1) follows by means 
of lemma 1 and lemma 6, where in this case we have 

Lemma 7: With ^f- = h , 

J^(a, a') = fc-T-Tfr + s)F 2 ( 1 + s; a, a'; 7, 7 - P\ \, j) (3.5) 

where F 2 ( 7 + s;a, a'; 7 ,7 - p; £, ^) is now computed by means of lemma 6 for s and 
p positive integers, and Re(h) > . Furthermore (1.2) follows directly for the special 
case s = p = . 

We may apply a similar proof of lemma 6 to prove the following result concerning 
the calculation of function F 2 . 

Lemma 8: For s and d positive integers and s>p, 

1. In terms of F 1 , we have 

F 2 {c + s;a,a';c,c + p; = (1 - ^ 'x 

g (p-s) (a') fcp F( s m + a ,. * W (3 ' 6) 

^ (c + p) m m\ k' k — h (h—k')(h — k) 

2. In terms of 2 Fi , we have 

F 2 (c + s;a, a>- c, c + p; * £) = (1 - \ )-(! - £)-' g - A)-™ x 

ft ft ft ft ^ (c + p) m m! fc' 

(3.7) 

E m (a) r (— s — m) r ,„ /i x _r r, / / kk! , 



r=0 



w/iere |fc| + \k\ < \h\ . 
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In particular, if s = p we have 

k V k V i k kk' 

F 2 (o + s, „, c,c + , -, -) = (1 - j)-(l - F,(a; -,,a';o; _. (3.8) 

Furthermore, from lemma 8, we have 

k k\ ( k\~ a ( k'\- a ' 
F 2 (c + s;a,a;c,c + s;-,-)=[l--j 

(3.9) 

E(a) r (—s) r / ^\~ r „ / , fcfe' . 

^H^H 1- *) A(tt+r;B;c+r; M)(M ) 



We can also apply Lemma 2 to obtain some new results concerning the calculations 
of the Appell's hypergeometric function F 2 . 

Lemma 9: For \k\ + \k'\ < \h\ , we have 

,-, / , k k 1 (a)r(c — a) r (d)2r ( kk'\ r k + k' ,„„„. 

Fa(d; °-° ;C - C; ^^ = g ( C ) r ( C )i r! (-7^) 2 ^ + 2r,a + r;c + 2r;^) (3.10) 

In particular, if k' = -k 

k k d d + 1 rr+lfc 2 

F 2 ((i; a, a; c, c; - , - -) = 4 F 3 («, c - a, - , — ; c, - , — ; ^) /or 2|fc| < (3.11) 

Proof: Form the following identity for the product of confluent hypergeometric func- 
tions [26] 

iFi(o; c; x) 1 F 1 (a; c; y) = ( ° )r(c ~ ° )r (-sy)^Fi(a + r; c + 2r; x + y) (3.12) 

we replace ir and y by ^ and ^ respectively, then integrate with respect to t 
both sides after multiplying through by t d ~ 1 e~ ht . The lemma follows immediately by 
applying (2.2) and (2.4) to the both sides. 

We can proof the particular case by re-placing k' with —k in (3.10) and applying 
the Gauss's duplication formula 

and then using the series representation expansion for 4 F 3 as given by (1.3). □ 
The lamme can be now used to show that: for Re(d) > and 2\k\ < \h\ , we have 

oo 

J t d - 1 e- ht 1 F 1 (a;b;kt) 1 F 1 (a;b;-kt) dt = ^ 4 F 3 (^ ^±1 a, b - a: b, \, ^) (3.13) 

o 
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or equivalently 

oo 

f (] i htm/ , b b+1 k 2 t 2 , , T(d) „ ,d d+ 1 , , b b+ 1 k 2 , . A , 

J ^ e - ht 2 F 3 {a,b-a;b, -,—;—) dt = -j^ 4 F 3 (-, —, a,b-a; b, -, —; ^) (3.14) 

o 

as a result of Ramanujan's theorem: 

b b + 1 x 2 

1 F 1 (a;b;x) 1 F 1 (a;b;-x) = 2 F 3 (a, b - a; b, -, — — ; — ). (3.15) 

The recurrence relations of the confluent hypergeometric functions can be used, 
by means of lemma 1, to develop several recurrence relations for the Appell's hyper- 
geometric function F 2 . For example, using the recurrence relation 

(b- a) 1 F 1 (a - l;b;z) + (2a - b + z) 1 F 1 (a; b; z) - aiF^a + l;b; z) = 0, (3.16) 

then, for Re(d) > -1 , the recurrence relation for F 2 follows 

£(d - l)F 2 (d; a, a'; b, b'; |, ^) = (a - b)F 2 (d - 1; a - 1, a'; b, b'; |, j) 

+ (b-2a)F 2 (d-l;a,a';b,b';^,^) (3.17) 

k k' 

+ aF 2 (d - 1; a + 1, a'; 6, b'; -, — ). 
In particular, if d = b' + s where s is positive integer, we have 

£(6' + a - l)F 2 (b' + s; a, a'; 6, 6'; ^, ^) = (a - 6)F 2 (6' +« - 1; a - 1, a'; 6, 6'; |, ^) 

+ (6 - 2a)F 2 (6' + a - 1; a, a'; 6, 6'; |, ^) (3.18) 
+ aF 2 (fc' + s - 1; a + 1, a'; 6, 6'; -, — ). 
Consequently, if a = 1 , the following recurrence relation follows 

^F 2 (6' + 1; a, a'; 6, 6'; \ k -) = (a b)F 2 (b>; a 1, a'; 6, 6'; |) 

+ (&-2o)F 2 (6';a,o , ;6,6';^) (3.19) 
fc fc' 

+ aF 2 (b';a + l,a';b, b';^,j)- 

and finally we have 

» m + !;«,«';«; *, *, _ („ - „) (1 - - £ )->,(„ - ,, „ - - ^ _ - ) 

+ (t - 2«)(1 - S )-(1 - ^)- rfl^^ti ,„„,„_,,. ) 



[h-k)(h-k')' 

k V , kk 1 



(3.20) 
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by means of (3.4). There are others recurrence relations for F 2 that can be develop 
in similar way, we mention three more: 

^±L^F 2{d . a , a >. b+h y +1 - * *1) = F 2 (d-1; a, a'; b, b'+l; * £)-F 2 (d-l; a-1, a'; 6, b'+l; J, £), 
/i c 6, ft /i /i n 6 

(3.21) 

(a + 1 - 6)F 2 (d; a, a'; b, 6'; ^) = aF 2 (d; a + 1, a'; 6, 6'; |, ^) - (b - l)F 2 (d; a, a'; b - 1, 6'; j). 

(3.22) 

fc fc' b(b - 1) fc ^ ^ , , A A' 6(6-1^ ,, , fc fc' 

F 2 (rf;a,a;6,6;-,-) = ^- r -F 2 (d-l;a,a;6,6;-,-)-^- r -F 2 (^ 

fc fc' 

+ {b- a)F 2 (d; a, a'; 6+1,6';-, — ) 

(3.23) 

4. Applications to perturbation theory 

Besides the important application of the integral (1.6) in computing the radial matrix 
elements of the electromagnetic interaction between the states of a relativistic electron 
in the Coulomb field, another important applications arises when we deal with singular 
Hamiltonians [27]- [30]. First we consider the following lemma. 

Lemma 10: 

1. F 2 (d; -n, -m; 6, 6'; x, 1) = ^^ 3 F 2 (-n, d, 1 - b' + d; 6, 1 - 6' + d - to; a:) 

2. F 2 (d; -n, -to; 6, 6'; 1, y) = £^ 3 F 2 (-m, d, 1 - b + d; b', 1 - b + d - n; y) 

3. F 2 (d; -n, -to; 6, 6'; 1, 1) = ^g^ 3 F 2 (-n, d, 1 - b' + d; b, 1 - 6' + d - to; 1) 

The proof of this lemma follows by means of the series representation (2.1) of the 
series F 2 , therefor we omit it. 

4.1. Singular Potentials 

A family of quantum Hamiltonians known as generalized spiked harmonic oscillators 
is given by the general Hamiltonian operator [27] -[30] 

H = H + ^, H = —^+x 2 + ^ (A>0,a>0,^>0) (4.1) 

acting in the Hilbert space L 2 (0, oo) . Eigenfunctions V e £ 2 (0, oo) of H satisfy the 
Schrodinger equation 

-?//' + <U 2 + 4 + — W = Eil> withV(0) = 0. 
^ x A x a ) 

The function ip is an eigenfunction corresponding to the eigenvalue E and the con- 
dition tp(0) = is called a Dirichlet boundary condition. It is known that the un- 
perturbed Hamiltonian H admits exact solutions by means of the Gol'dman and 
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Krivchenkov wavefunctions 



-n; 7 ;x 2 ) for n = 0,1,2,.... (4.2) 

where 7 = 1 + \\/l + AA, ; and the exact eigenvalues are given by 

E n = 2(2n + 1 ), rc = 0,l,2,... (4.3) 

The orthonormality of the Gol'dman and Krivchenkov wavefunctions {ip n (x)}%L as 
well the computation of the matrix elements of the perturbed operator V(x) = x~ a in 
terms of {ip n (x)}%L follows immediately by means of lemma 10. Indeed, we have 

/>oo 

/ ip n {x)ip m (x) dx = 6 mni n,m = 0,1,2,... 
Jo 

as a result of 

OO 

/ x^- 1 e- hx \F 1 (~n;r,hx 2 ) 1 F 1 (-m;r,hx 2 ) dx = (4.4) 
J 2 /i 7 (7)„ 



The proof of the last integral reduces, after simple change of variables, to the com- 
putation of 

F 2 (r,-n,-m;'Y,r,lA) = - 2F 1 (-n,j + k; T , 1) = —i— V ^ "^J ^" = -r^-5 mn , 

fc=0 l7j ™ fe=0 l7jn 

since the product (-m) fe (-fc)„ differs from zero only for n = k = m . Furthermore, the 
matrix elements of the singular operator V(x) = x~ a follow immediately by use of 
Lemma 10, and this yields 

Lemma 11: // 27 > a , then for all pairs of non-negative integers and m and n we 
have that 

00 

f ^- a ~ 1 e-^iF 1 {-n^,h^) 1 F 1 {-m^,hx l )das = 

(4.5) 



^- 7 (f) n r(7-f) 



a _ a a 



2 ( 7 )„ M- mn --A--^l---n-l). 

Consequently, the matrix elements assume the explicit forms 



which can be used directly for a variational study of the perturbed Hamiltonian H . 
Note, in particular, for n^O 



<^o|x- Q |V„ >= (-1) 1 



(7)» r(7-f)(f)n 
n\ r( 7 ) ( 7 )„ ■ 
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4.2. Perturbed Kratzer Potential 

There has been great interest in the Coulomb potential with polynomial perturba- 
tion from both mathematical and physical points of view [31]. However, the Kratzer 
potential [32] 

V(r) = ~ + ^, (B>0,A>0) (4.7) 

which contains the Coulomb potential as a particular case (A = 0) has apparently not 
been studied in such a context [31]. The radial part of the Schrodinger equation for 
Kratzer potential reads 

d 2 2d 1(1 + 1) ( B A\] . . . . 

* V>(r) = Eip(r) (4.8) 



dr 2 r dr r 2 \ r r 2 

The (un-normalized) wavefunctions are given explicitly by 

Mr) =r fl e-3^TT r 1 F 1 (-n,2a + 2,— ^—r), n = 0,1,2,... (4.9) 

n + s + 1 

where 

s = - l - + \^AA+{2l + l) 2 . (4.10) 
The exact eigenvalues are given by the expression 

E m = - T( n2 , n,/ = 0,1,2,... (4.11) 

4(n + s + 1)^ 

For the case of the Coulomb potential A = and s = Z, where I = 0,1,2, ... is the 
angular momentum quantum number. The normalization constant (with the angular 
factor 4-7T omitted) appropriate for the wavefunction solutions (4.9) can be computed 
by means of the condition 

/>oo 

/ \ip n i\ 2 r 2 dr, 
Jo 

which leads to 

C-J = ^ r 2 ^ 2 e-^[ 1 F 1 (-n,2 S + 2,—^—r)] 2 dr 
Jo «tst J 

The integral just mentioned can be computed by an application of Lemma 10, which 
yields 

C- 2 = 2B- 2 °-*{s + n + l)^rg£ + g n , (4.12) 

If the potential (4.7) is perturbed by the operator Xr a , a > , that is to say if we 
consider 

B A 

V(r) = -- + — + \r a , A>0,a>0 

The matrix elements for the operator r" by means of the (normalized) wavefunctions 
(4.9) and (4.12) are given through the integral 



/>oo 

< <Mr Q |Vw > = C nl C ml / r 2s+2+a e-^^TT + ^TT> 

Jo 



X 



Br Br 

iFi(-n; 2s + 2; — _) 1 F 1 (-m; 2s + 2; ■ — )dr 

n+s+1 m+s+1 



Integrals containing confluent hypergeometric functions . . . page 15 



which yields by means of the series representation (2.1) of F 2 that 

-a-2s-3 



< Vw|r Q |Vw >= T(2s + 3 + a) C nl C ml 
(2s + 3 + a) k (-n) k ( 2m + 2s + 2 



D 



+ 



fe=0 



(2s + 2) fe k\ \2s + (n + m) + 2 



2 \n + s + l m + s + 1 

k 2n + 2s + 2 
2^1 (2s + 3 + a + k, -to; 2s + 2; -). 



2s + (n + m) + 2 ' 
(4.13) 



In particular, the diagonal matrix elements (i.e. m = n) are given by 

< *,n*. > - i*~*<» + . + ■r 2 - r( ^::y )( - i n r ) " > 

3 F 2 (-n, 2s + 3 + a, 2 + a; 2s + 2, 2 + a - n; 1). 



(4.14) 



The matrix elements (4.13) and (4.14) can also be used immediately, for example, for 
a variational study of the spectrum of H = -A + V(r), where V(r) is a power sgn(q)r a . 



5. Conclusion: 



We have shown that a simple integral can be used to generate many identities that are 
very useful in mathematical and physical applications. It is clear that the confluent 
hypergeometric function iFi(-n,a + l,x) can be expressed in terms of the associated 
Laguerre polynomials L%(x) through the identity 

1 F 1 (-n;b+l;x) = -^ T ^-L b n (y) 

and many of results dealing with the applications can be expressed equally well in 
terms of L%(x) . Some interesting new identities for class of intergrals containing pair 
of associated Laguerre polynomials L%(x) with weight measure {t d ~ 1 e~ ht } can easily 
be derived on the basis of the explicit results presented here. 
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Appendix I 



In this appendix we collect some results that follow by means of integral (2.4), 
and may be useful for further applications. Indeed, the possibilities implied by inte- 
gral (2.4) are not limited to the present results; many other identities follow by use 
of the available expansions of the Gauss hypergeometric functions 2 Fi in terms of 
elementary functions; see for instance [33]. 
For Re(d) > and \h\ > \k\ 

oo _ 

J f-V'xiMM;**) dt= M 2Fl (d,a;d;l) = M lFo (a; ^ = £^ ^ _ (L i) 
o 

For Re(d) > , Re(h) > and a + d + 1 ^ 0, -2, -4, . . . 

F (a- a + d + l . ht ]Hf - r W p (n d . a + d + i . u r(d) r(i)r(^f±i) 



For i?e(a) > , iie(ft) > and a + b ± -2, -4, -6, . 



2 '2 ; ^ d r(i±^)r(i±^)' 



= 2 



r(a) r(i)r(^ + i) 



(1.3) 



a — b 



r(f)r(^±i) r(|)r(s±i) 



For i?e(2a) > and > |fc| , 

OO 

J t^^e-^xF^a + ^a-^t) dt = ^^ 2 Fx{2a, a + 1; a; |) = r(2a) (/i + fc)(/i - fc)" 20 " 1 . (1.4) 
o 

For i?e(a) > -1 and > |fc| , 



OU 

/ *°e-^ 1 F 1 (2a; a; fct) dt = T ±±ll (l + £) (l - £ 
o 

For i?e(a) > and \h\ > \k\ , 

OO 

o 

For Re(a) > -\ and \h\ > \k\ , 



-2a-l 



(1.5) 



; ^(a+^a;**) dt = ^ 2 F 1 (a,a+l;2a;|) = M 





1 

2 


'1 


I— 1 




['"5 








A-s] 






2~ 


h 2\ 





OO 



e _/lt iFi(a;2a;fct) dt 



r(a + ^) 



1 1 

2 + 2 



l-2o 



l-2a 



(1.6) 



(1.7) 
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For Re(a) > and \h\ > \k\ , 



J t-^-^iF^o-^ajfct) d*=^r 2^1 (a, a 2a; |) 



11/ fc' 



l-2a 



(1.8) 



For Re(a) > \ and \h\ > \k\ , 



oo 



3e"' lt iFi(a;2a;fct) di 



r( a -|) 



1 1 / fc 
2+W 1 ^ 



l-2a 



(1.9) 



For i?e(a) > and \h\ > \k\ , 

OO 

o 

For Re(a) > -\ and > |fc| , 



oo 

/«.-». 



e~ M 1 F 1 (a;2a + l;kt) dt = 



r(q + |) 

h a +i 



1 1 / k 
2+W 1 ^ 



2 2 V 



-2a 



For i?e(a) > -1 and \h\ > \k\ , 



r(a + 1) „ , _ 1 _ fc, 



y t°e- fct iFi(a+^;2a + l;fct) = - 2 Fi(a + 1, a + -; 2a + 1; ^) 



r(a + i) 



1 - 



1 1 / k 

2 + 2V-h 



For Re(a) > -\ and > |fc| , 



^e~ ht - L F 1 (a + 1:2a + l;kt) dt = 



h a +i 



1 1 / k 
2 + 2i l -h 



For Re(a) > -\ and |A| > , 



oo 

J t a -h- ht 1 F 1 (a;2a + l:kt) dt = 



r(q+|) 

h a +i 



1 1 / k 

2+W 1 ^ 



-2a 



For i?e(a) > and \h\ > \k\ , 



OO 



e _/lt iFi(a+ ^;2a+ l;fci) dt 



1 



£(a) 



1 ll k 

2 + 2r--h 



-2a 



For \h\ > \k\ and J?e(a) > -1 



t-e-^a^fci) dt=%±il 2 F 1 (a + 1,1; 2; J) - r( "> 



/i Q (h-k) a 



-i -2a 



(1.10) 



(1.11) 



(1.12) 



(1.13) 



(1.14) 



(1.15) 



(1.16) 
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and for \h\ > \k\ , 



oo 

/■ 



1 



r -ilog(l-|), if a = 



e- M l F l {a+l;2-kt) dt = - 2 Fi(l, a + 1; 2; -) = { 



h 



a ' a(h-k) a 



, if a ^ 



(1.17) 



In particular, for \h\ > \k\ , 

oc 

J e- /lt 1J F 1 1 (2;3;fct) dt = /i" 1 2 F 1 (1, 2; 3; |) = ~ 



(1.18) 



For the particular calculations of the hypergeometric functions 2*1 (a + 1, 1; 2; x) in 
terms of elementary functions, we refer to Ref. [15]. 

Direct aplications of (2.2) also yields: For Re(a) > and s > ^f- , we have 

o 

(1.19) 

In particular 



oo 

1 1 2 "-^ 



)Fi(-;o + 



1 fc 2 t 2 T 2 



2' 16 



eft 



r(2 Q ; 

Q 2a 



k 2 

2 Fi(a,a;2a; -~). 



(1.20) 



as a result of the identity 



1 z 2 

iFi(a;2a;z) = e* Fi(-; a + -; — ), 



which also leads to the particular case 



o*! (- 



3 k 2 t\ 



2' 16 



1 k 2 

^ = -^2 l0 S(l-^2)- 



(1.21) 



